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^•f-^ I We discuss the impact of nuclear effects on the higher twist terms using a particular 

example of the structure function xF^ as extracted from a QCD fit to neutrino deep 
inelastic data on a heavy nucleus target. 

^ ■ The operator product expansion is a common theoretical framework in analyses of 

^ . deep inelastic scattering (DIS) in QCD. The operators can be ordered according to their 

I twist yielding the series in for physical observables. For example, for the structure 

O I function F3, this reads 

9^'- xFsix,Q') = xFt^ix,Q') + ^^^^^ + ... (1) 

r-| ! The first term in this expansion (the leading twist, LT) dominates at sufficiently large 

momentum transfer and invariant mass W"^ = + Q'^{1 — x)/x. The LT structure 



X 



functions are constructed in terms of parton distribution functions (PDFs), which are 
universal for charged lepton and neutrino scattering and have clear probabilistic inter- 
pretation. An accurate knowledge of these plays a key role in the extraction of possible 
contributions of new physics at new collider energies, non-accelerator physics (cosmic neu- 
trinos) and, as observed more recently, in the interpretation of forthcoming high precision 
experiments on neutrino oscillation. 

PDFs are not directly observable and extracted from data in global QCD fits (see, e.g., 
[0]). The higher twist (HT) terms are presently poorly known and currently is a subject 
of both theoretical and phenomenological studies 0, |^, |I], ^, ^ . A better understanding of 
HT terms, in particular their role in describing low and high x DIS data, is important 
and provides valuable information on quark-gluon correlations inside the nucleon. 

Phenomenological studies of PDF and HT terms are affected by a number of uncertain- 
ties. These include the effects due to higher order perturbative terms in renormalization 
group analysis, threshold resummation effects, target mass corrections, the sensitivity of 
QCD fits to the choice of parameterization of PDFs and HT, etc. If nuclear data are 
involved into an analysis, nuclear effects will also have impact on the extracted PDF and 
HT terms. 

In the present contribution we discuss the size of nuclear effects on a particular example 
of the HT term in the structure function F3 as extracted from neutrino data of CCFR 



collaboration on iron target (for more detailed discussion see 0). In our fit, we separate 
the nucleon structure function F^{x,Q'^) into a sum of the LT and the HT terms, (|l|), 
and parametrize xF^"^ at some scale = Qq in terms of a simple function, 

xF^'^{x,Ql) = aix"^(l - x^'il + a^x). (2) 

Then we apply the renormalization group equation in order to calculate evolution of xF^^ 
with Q^. We solve the renormalization group equation in the leading (LO), next-to-leading 
(NLO) and next-to-next-to-leading (NNLO) logarithm approximations of QCD. In doing 
so we expand the leading twist structure function xF^'^ in terms of its Mellin moments 
within the framework of the Jacobi polynomial method and then apply the evolution 
equations to the moments (for more detail see |^ and references therein). 

We fit 116 data points on the structure function F^ from the CCFR experiment in 
the kinematical range of 0.0075 < x < 0.75 and between 1.3 and 200 GeV^. The fit 
parameters are 02, 03, and 04 of Eq.(0) at the scale Qq, the values of the function H^, at the 
center of each Xj-bin of the CCFR data set, as well as the QCD scale parameter The 
parameter ai was fixed by normalizing (|^) to the Gross-Llewellyn-Smith sum rule, which 
was calculated in QCD to the second order in as 0, 5'gls = 3(1 — as/n — 3.25(a5/7r)^). 

It should be noticed that in general the higher twist terms can be of two kinds: those 
which have the kinematical nature, e.g. the terms due to finite target mass, and those 
which arise due to higher twist operators and reflect the quark-gluon interaction effects 
in the target. In order to ensure that the HT term in (|I]) describes effects due to quark- 
gluon interaction in the target we explicitly take into account target mass correction by 
substituting the Mellin moments by the Nachtmann moments in the Jacobi polynomial 
expansion of xF^^ . 

In order to apply corrections for nuclear effects in our analysis we first calculate the 
"EMC ratio" for the iron target, R^{x,Q'^) = ^F^^x.Q"^) / F^ {x,Q'^), with F^" the 
structure function of the ^^Fe26 nucleus and F^ = ^(-Ff -|- F-^) the structure function 
of an isolated isoscalar nucleon. Then we extract "data" on the structure function of 
the isoscalar nucleon from the CCFR data as Fi^{x,Q'^) = F^'^^^{x,Q'^)/R3{x,Q'^). In 
calculating the ratio -R3 we observe that the bulk of neutrino data with > 1 GeV^ is 
located in the region of x > 0.1. For this kinematical regime it is usually assumed that 
nuclear DIS of leptons from nuclear targets can be viewed as incoherent scattering from 
bound nucleons. Major nuclear effects found in this region are due to nuclear binding 
and Fermi motion |^. The relation between a heavy nucleus structure function and the 
proton and neutron structure functions can be written as follows 

where x' = Q'^/2p-q is the Bjorken variable of the bound nucleon with the four-momentum 
p, q is four momentum transfer and 7 = \q\/qo is the 'velocity' of the virtual boson in 
the target rest frame. The averaging is done with respect to the nuclear spectral function 
and the last term in Eq.(|^) takes into account that the neutron, N, and the proton, 
Z, numbers are unequal in heavy nuclei. Note also that bound proton and neutron are 
off-mass-shell and their structure functions depend on the nucleon off-shellness p'^ as an 
additional variable that causes an additional correction . The nuclear spectral function 
and other details of the present approach are discussed in . 



Our results are shown in Fig.[T|. One can observe from Fig.|l| that the scale which 
determines the magnitude of the HT term is 1 GeV^. One also sees that the magnitude 
of the HT term depends on the level to which the perturbation theory analysis of F^^ 
is done. The more perturbative corrections are included into the evolution equation, the 
less room is left for the function H^lx). The separation of nuclear effects from data leads 
to a further suppression of the HT term H^{x). The effect of nuclear corrections is most 
pronounced at large x, where we observe a systematic reduction of i?3(a;) as compared 
with no-nuclear-effects analysis. Nuclear corrections result in the decrease of the values 
of the function H^i^x) at x > 0.6. We also found that the separation of nuclear effects 
causes the to decrease for about 10% leading to the shift of the value of as{Mz) for 
about 2 ■ 10^^. As a final remark, we comment that a reduction of the HT terms because 
of nuclear effects at large x has recently been observed in the analysis of charged lepton 
DIS data on proton and deuteron targets 0] .f| 
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" However, the HT term in the structure function F2 of ^ appears to be an order of magnitude less 
than the central points of H3 of Furthermore, the effect of strong correlation between the order of 

pQCD analysis and the magnitude of the HT term, first observed in y] for neutrino data, is not seen in 
the recent analysis of 
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Figure 1: The function H^{x) in the units of GeV^, which describes the strength of the higher twist 
term in the xF^ structure function as extracted from the fit to the CCFR neutrino data. The labels on 
the figure indicate the level to which the perturbation theory analysis of xF^'^ is done. Only statistical 

errors are shown. 



